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We present a formula that allows one to calculate the pion form factor in the timelike region 2m 7V < 
t/s < 4m T in lattice QCD. The form factor quantifies the contribution of two-pion states to the 
vacuum polarization. It must be known very accurately in order to reduce the theoretical uncertainty 
on the anomalous magnetic moment of the muon. At the same time, the formula constitutes a rare 
example where, in a restricted kinematic regime, the spectral function of a conserved current can 
be determined from Euclidean observables without an explicit analytic continuation. 
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The ab initio calculation of hadron properties in lat- 
tice QCD is a mature and ongoing research activity. For 
a compilation of results for some quantities of high phe- 
nomenological impact see Ref. One important as- 
pect of lattice QCD is that it is formulated in Euclidean 
space. Therefore the form factors are necessarily ex- 
tracted for momentum transfers in the spacelike region, 



q 2 = q 2 — q 2 = —Q 2 < 0. At first sight, it thus seems that 



a form factor in the timelike region, or more generally any 
quantity which involves the notion of 'real time', is only 
accessible via analytic continuation. The latter typically 
involves solving a Fredholm equation of the first kind 
such as Eq. ((3]), which represents a numerically ill-posed 
problem 0. 

The goal of this letter is twofold. First, to provide an 
example where, in a certain kinematic regime, a spectral 
function (to be defined shortly) can be extracted from 
the Euclidean field theory without an explicit analytic 
continuation. And secondly, to lay the ground for future 
calculations of the pion form factor in the timelike re- 
gion, which may impact the theoretical determination of 
the anomalous magnetic moment of the muon. Given the 
high likelihood of a new (g — 2) M experiment at Fermi- 
lab Q , this application is particularly timely. 

For the electromagnetic current j° m = |w7 M u 



■>H 3 ' 

the spectral function is defined as 



Av(fc) = — 

Z7T 
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Due to current conservation and Lorentz invariance, the 
tensor structure of p^ v is 



p(k 2 



(2) 



The spectral density p is non-negative. The lead- 
ing hadronic contribution to the vacuum polarization 
e 2 H(Q 2 ) in the spacelike domain can be expressed 
through it via a once-subtracted dispersion relation, 
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The function H(Q 2 ) can be calculated in lattice QCD 0- 
7]. On the other hand, p(s) is related to experimental 



g(f) 
12tt 2 ' 



a(e+e -> hadrons) 
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The denominator is the treelevel cross-section e + e~ — > 
p + p~ in the limit s 3> m 2 ^ , and we have neglected QED 
corrections. At low energies, the spectral density is given 
by the pion form factor defined in (|12|) Q, 



P(s) 
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|^(0)| = 1. (5) 



The relation holds for 2m,r < y/s < 3m w , and even up 
to 4m w if the electromagnetic current is replaced by the 
isospin current in the definition of p(s). 

In infinite volume, the spectral density is a continuous 
function above the two-particle threshold. In finite vol- 
ume, where simulations are carried out, it is a collection 
of delta functions. What is needed is an explicit formula 
that relates the energy and matrix element of one indi- 
vidual state living on the 3d torus to the infinite-volume 
spectral density at the same energy. In this letter we 
derive the formula 



IF^E)] 2 = (q<f>'(q) + k 
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Here E equals the invariant mass of the two pions, fc is 
related to E via Eq. (fT3")) . Si is the scattering phase shift 
in the unit isospin, p-wave channel and is a vector- 
current matrix element between the vacuum and a two- 
pion state of energy E on the torus, see Eq. (p2 
Finally, q = |p and (f> is a known kinematic function 
The scattering phase 5\{k) can be extracted (see [10l. 
and Refs. therein) from the finite- volume spectrum using 
the Liischer formula HG3, E< 1- 13) • 

The derivation of Eq. ^ is closely related to work 
by Lellouch-Luscher [13[ on the matrix element deter- 
mining the K — > 7T7T decay rate. Eq. ([SHS]) show that 
this formula, remarkably, allows one to extract a spectral 
function from Euclidean field theory without an explicit 
analytic continuation. 
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II. PRELIMINARIES 

We will assume here that isospin is an exact symmetry 
of QCD, and focus on the vector isovector channel 

7=1, J PC = \ — . (7) 

Both in experiment and in the Euclidean theory, the sym- 
metry of the final state can be selected. Therefore, in 
spite of the physical photon not having definite isospin 
quantum numbers, we will consider the case of a gauge 
boson coupling to the isospin current 

jZ = jrf^, ^=(d)- (8) 

In infinite volume, and with quark masses set to their 
physical values, the states in the symmetry channel ([7]) 
are necessarily two-pion states for center-of-mass energies 

Im-x < E < 4m n . (9) 

We are thus led to consider pion scattering in this energy 
range. First, the normalization of single-pion states is 

(k a\k' b) = 5 ab {2Tr) 3 2E k 5(k - k'). (10) 

The Tj matrix describing the scattering of two pions is 
defined in full generality in (Q, Sec. 2.1). It carries 
an isospin index / = 0,1,2; here we are interested in 
the 1=1 channel. The T\ matrix can be expanded 
in partial waves with amplitudes tjg , with corresponding 
phase shifts Sn . We restrict ourselves to the I — 1 partial 
wave in the elastic regime, 

T\ = 487r£Wtii costf, t u = -. (11) 

From now on the phase shift will simply be written Si. 
When the pions couple to the photon, they can be pro- 
duced for example in the reaction e + e~ — > 7* — > rnr. 
With |(vrp7Tp') a ) = ^^•| 7r p 7r p')) the pion form factor in 
the timelike region is defined as [l4| 

(0|i°|(7r p 7v) 6 ,in) = -<(7r p7 r P b ,out|/|0) (12) 
= 5 ab V2i{p' -p)F w (E wn ). 

The two-pion states on a three-dimensional torus of 
dimensions L x L x L at vanishing total momentum have 
been studied in detail in Q. The vector states (i = 1) 
are found exclusively in the Tj irreducible representation 
of the cubic group. Their norm will be taken to be unity. 
We denote by J5 W the energy of one such state. The 
effective momentum of the pion k v is then defined by the 
equation 

E m = 2^ml + kl. (13) 

The discrete values of fc^ = k in the box and the 
infinite- volume phase shifts are related by 0, [HJ 

kL 

Si(k) + (t>(q) =n7T, 9= IT' ( 14 ) 



The function <f>(q) is defined by tan</>(g) = — ^(uq'i) > 
where Zoo(l;q 2 ) is the analytic continuation in s of 
-2oo(s;<? ) = ^E n6 z 3 ( Tl 2_ ? 2 )s ■ 



III. QCD COUPLED TO SU(2)] GAUGE 
BOSONS IN THE BROKEN PHASE 



The theory we consider is QCD with at least two de- 
generate light flavors of quarks u and d. The other flavors 
are assumed to be sufficiently massive that the only unit- 
isospin states in the interval (|9j) are the two-pion states. 
This is realized for physical values of the quark masses. 

We now wish to establish Eq. (jSJ) in the regime where 
the gaps between finite-volume energy eigenstates are 
substantial, which is the situation that can realistically 
be achieved in Monte-Carlo simulations. The idea is to 
couple the quarks infinitesimally (the coupling will be de- 
noted by e) to an external vector particle with mass in 
the range ((9]). On one hand, when it is degenerate with a 
two-pion energy eigenstate in the box, a level splitting of 
order e occurs; on the other hand, in infinite volume the 
resonant production of the vector particle in pion scat- 
tering leads to an 0(e) change in the phase shift. The 
phase shift and the energy levels in the box must be in 
correspondence through the Luscher formula flT^]) both 
before and after switching on the coupling e. Eq. ([6]) fol- 
lows from the difference of these relations. The proof is 
thus almost identical to the proof of the Lellouch-Liischer 
formula for the kaon decay |13| . The only qualitative dif- 
ferences are that we are in a different symmetry channel 
and, more importantly, that we invoke an external par- 
ticle. A similar formula for a matrix element involving a 
two-particle state was derived in [l5j , in the large- volume 
regime where individual states are too narrowly spaced 
to be individually resolved. 

A concrete model for the massive vector boson is not 
required in the proof, but it may be reassuring for the va- 
lidity of the argument that the situation sketched above 
can be realized within a renormalizable field theory. A 
specific realization, then, involves gauging the isospin 
symmetry (this was, incidentally, the original motivation 
of Yang and Mills [16[ to introduce non-Abelian gauge 
theories) . There are then three gauge bosons of the group 
SU(2)i. Since we want them to be massive, we assume 
that the gauge group is spontaneously broken by a Higgs 
mechanism at a scale well above Aqcd, the scalar field 
being in the fundamental representation of SU(2)j. The 
three gauge bosons W a then form a degenerate triplet of 
mass M, which we assume to be in the range They 
have the same quantum numbers as the p meson and 
couple to the isospin current. 



3 



A. Two-pion energy levels in finite volume 

We start with the system in finite volume and study 
the effect of the massive gauge bosons on the low-lying 
spectrum, which consists of two-pion states. The gauge 
boson field operator reads 



AUx) 
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), (15) 



with 



,] = S ab S aa 'Skk' and the polarization vec- 



tors eZ are such that at k 



0, e° n 



and e° 



Q — V OllU — 

The coupling to the massive gauge bosons is treated as a 
perturbation, with 



H int = +e J j a »{x)Al(x). 



(16) 



The box size L is now chosen such that the energy 
eigenvalue of a unit-norm, non-degenerate state °f 
the two-pion system overall at rest in the box satisfies 
= M when e is set to zero. Upon switching on 
the SU(2)j gauge interaction, we apply degenerate per- 
turbation theory of quantum mechanics. The two-pion 
state mixes maximally with the massive gauge boson 
state \W%) = Oq^IO), since they are initially degener- 
ate, and the splitting between the two resulting energy 
levels is given by the off-diagonal matrix element, 



E% r = M±\A\ 



(17) 



where q 2 - M 2 = ±2M|.A| + 0{A 2 ). Using the definition 
of the pion form factor in the timelike region (|12j) and the 
phase F n — l-F^e 4 * 1 due to Watson's theorem, we obtain 



2 sab 



ATi 



e 2 6 



M\A\ 



\F„(E„„ = M)\ 2 4k 2 cos6e 



(22) 



Next we translate this into a change in the phase shift 
via Eq. (TIT]) . 



A5 1 = T 



L3 

2 7T 



12irM\A\' 1 E, 



(23) 



C. Connection between the energy levels in 
finite-volume and the scattering phase shift 

We now go back to the finite volume system where L is 
tuned so that E n7r = M. Upon switching on the pertur- 
bation (|16p . this degenerate energy level in the box splits, 
corresponding to k — > k + Afe, where in view of Eq. (|13f 
[T7|) , Ak = ±A ^F- . The phase shift that corresponds 
to the perturbed state by Liischer's formula (TlH) differs 
from the phase shift corresponding to the unperturbed 
state by 

S 1 {k)^S 1 (k) + ^^-Ak + AS 1 (k). (24) 

In Eq. (1141) . this variation must be compensated by the 
change in the second term of the LHS, 



</>(<?) -> <P(q) + <f>'(q)Aq, Aq = q 



Ak 



(25) 



with 



mH int \W2,) =8 ab 8 aa ,A 



(18) 
(19) 



Altogether, we then obtain the condition 



A8 1 (k) = -(q ( f>'(q) + k 



dSi(k)\ Ak 
dk 



(26) 



Inserting expression ([2"B"| for the change in the phase shift 
due to the resonant scattering, one then obtains Eq. ([5]). 



B. Pion-pion scattering in infinite volume 

The coupling of quarks to the massive gauge bosons 
affects their scattering amplitude. At ordinary energies, 
the effect is 0(e 2 ). However, at the energy (TT7)) the effect 
is enhanced to 0(e), due to the resonant production of a 
massive gauge boson, whose propagator reads 



iS ab 



M 2 



ie V M 2 



(20) 



The change of the scattering amplitude due to this pro- 
cess then reads, at the energy (fTT|) and with center-of- 
mass frame kinematics, 



(21) 



AT\ ^(7Tp7r_p) a -> W^(Tr p ,n- p , y 
= ((^vr_ P O b ,out|j c |0) 



q 2 - M 2 



i|j d |(7r p 7r_ p ) a ,in), 



D. Consistency check in the absence of pion-pion 
interactions 

In the case where the two pions do not interact, the 
71 th two-pion energy level E^ passes through M for [l3[ 



2tt 

L = — V n : 



1 < n < 6, 



(27) 



kn being related to E n7r by Eq. (Ti~3"]) . From the definition 
of 0, one sees that for q 2 — > n an integer, 



b'(q) = (27T) 



2<t 



(28) 



where v n is the number of vectors z £ 1? with z 2 = n. 
The result © then becomes for non-interacting pions 



\ F ^ = M)\* = ^^\A^. 



(29) 
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Finite-volume matrix element 
5 11: Eq.(32) — 



IV. ILLUSTRATION 



§11=0 




FIG. 1. Finite- volume matrix element \A^,\ 2 calculated (a) 
using Eq. ([6jl and (b) using Eq. (|29[l . where the pion inter- 
actions are neglected. The phenomenological input for 5i(k) 
and F 7T {E) is given in Eq. ([32]) and ([33j). The values of the 
ground state energy is indicated for three different box sizes. 



One can check this last equation directly. Take for 
instance as a two-pion state on the torus \ip%) given by 



Here we apply Eq. ^ to infer the dependence of the 
matrix clement A^, on L from experimental scattering 
data. The scattering phase 5\ (k) is parametrized by [§] 



cot Si(k) 



4*£ 



i-2 

P 



(32) 



with k n 



4ml. 



t, p — 2 y "V ~~ ■"'v I n ^ s effective range formula, 
the scattering phase is thus determined by the mass and 
width of the p meson. Secondly, we extract the pion form 
factor from experimental data compiled by the PDG [l?} ■ 
At low energies the data is well described by 

\FAE^)\ 2 = vo + v 1 ([(E^/m„) 2 -v 2 ] 2 + v 2 3 y\ (33) 
where the fit parameters for the interval 2.Qm n < E n „ < 
4.4m T take the values vo = 0.6473, V\ — 10.59, v$ — 
0.1271, and the x 2 /d.o.f. amounts to 1.2. We did not fit 
V2-i but set it to the value m p jm^ — 5.553. 

Using formula Eq. (j6]), we can predict the magnitude of 
the finite-volume matrix element A^p up to the inelastic 
threshold in the 1=1 channel (4m 3r ). The result of this 
exercise is displayed in Fig. (JTJ, along with the finite- 
volume matrix element one would expect if there were no 
interactions between pions. As already discussed in Q, 
the size of the matrix element \A^, \ 2 is of order L~ 3 in the 
absence of interactions between the pions, but it would 
be strongly enhanced in the vicinity of a resonance. 



CONCLUSION 



,abc 



m = —a b Ja%\0), 



S ab S r 



(30) 



p<=3 and a°' a pion creation operator. 



Thus from 



withp - ... ., ,. p 

the difference in normalization of the finite and infinite- 
volume states (Eq. [3D] and ITU]) one predicts 



8kl\F n \ 2 = M 2 L 3 \A^ 



(31) 



which matches ([29]) for v n = 6. This is precisely the 
multiplicity v n of n = 1 and 4 for which the only vectors z 
are (y/n, 0,0) and images thereof under cubic symmetries. 



We have derived a formula which connects Euclidean 
observables with the vector current spectral function in 
the elastic regime. Such a relation could have interest- 
ing applications in other systems, such as cold Fermi 
gases [lg|. Also, the idea to invoke an external particle 
weakly coupled to QCD is likely to lead to other impor- 
tant relations. Apart from its theoretical appeal, Eq. ^ 
provides a more direct way to compare lattice calcula- 
tions with experimental determinations of p(s) than the 
dispersion relation ([3]). To complement the low-energy 
experiments constraining the leading hadronic contribu- 
tion to (g — 2) M , lattice simulations would have to reach 
very low pion momenta. For this purpose [l9l - l2l| , as well 
as to obtain the factor dSi(k)/dk [221, twisted boundary 
conditions could prove very useful. 

I thank H. Wittig and D. Berneckcr for helpful discus- 
sions and M. Liischer for comments on the manuscript. 
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